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Abstract 

Let Bi, B2 be balls in finite-dimensional real vector spaces E\,E2, centered around unit length 
. vectors v\ , V2 and not containing zero. An element in the tensor product space E\ (g> E2 is called 

Bi ® inseparable if it is contained in the convex conic hull of elements of the form wi ®W2, where 
wi 6 Bi, W2 € B2. We study the cone formed by the separable elements in Ei®E2- We determine 
■ the largest faces of this cone via a description of the extreme rays of the dual cone, i.e. the cone of 

5^ ' the corresponding positive linear maps. We compute the radius of the largest ball centered around 

vi ®V2 that consists of separable elements. As an application we obtain lower bounds on the radius 
of the largest ball of separable unnormalized states around the identity matrix for a multi-qubit 
system. These bounds are approximately 12% better than the best previously known. Our results 
are extendible to the case where Bi, B2 are sòlid ellipsoids. 

(N . 

1 Introduction 

, Let K, K' be closed convex pointed cones with non-empty interior, residing in finite-dimensional real 
t— ( ' vector spaces E, E' . Then an element w £ E <g) E' of the tensor product space is called separable if it 
can be represented as a convex combination of product elements v®v', where v S K , v' E K' . It is not 
hard to show that the set of separable elements is itself a closed convex pointed cone with non-empty 
interior. This cone is called the K ® íí'-separable cone. Cones of elements that are separable with 
respect to more than two initial cones are defined similarly. 

Separable cones have many applications in Mathematical Programming. So, the dual cones to the 
cones of positive polynomials, which frequently appear in optimization problems 0,0 , can be repre- 
sented as separable cones. In Quantum Information Theory the set of unnormalized unentangled mixed 
states of a multi-partite quantum system also forms a separable cone. In this case the underlying cones 
are cones of positive semidefinite matrices. These separable cones and the corresponding positive maps 
have been subject of intense study in the recent Quantum Information Theory literature |5] . [5] . [TU] . but 
drew attention of the mathematical community also before the emergence of this field |T|.|12 | .[T3 | .|16 | . 

In this paper we treat ball-ball separable cones, i.e. cones of K (g> if'-separable elements where the 
underlying cones K, K' are conic hulls of closed Euclidean balls or sòlid ellipsoids not containing the 
origin. Such cones have a relatively simple structure. Thus they are suitable for the approximation of 
more complex separable cones. This can be done by approximating the underlying cones by appropriate 
ball-generated cones. The idea of replacing an underlying cone by a ball-generated cone was put forward 
by Leonid Gurvits and Howard Barnum, who successfully used it to obtain lower bounds on the largest 
ball of unnormalized separable elements around the identity matrix for multipartite systems 0,0. In 
this contribution we compute several characteristics of ball-ball separable cones exactly, which allows 
for a more cfficient application of such approximations. 

One such application makes use of the fact that the cone of positive semidefinite hermitian 2x2 
matrices is isomorphic to the 4-dimensional Lorentz conc £4 and hcncc is also a ball-generated cone. 
This enables us to refine Gurvits' bounds for the case of multi-qubit systems. We prové that for a 
3-qubit system a ball of radius \J 4/5 around the identity matrix consists only of separable elements, 
as opposed to the best bound -y/8/11 known previously [Hj- For systems consisting of more than 3 
qubits we obtain an improvement of roughly 12% with respect to the best bounds known before 6 . 
Namely, we prové that for an m-qubit system, a ball of radius —j = around the identity matrix 

consists only of separable elements. The exponent in the asymptotics (as m — > 00) of this bound is 
the same as the one obtained in . Recently Stanislaw Szarek showed that this exponent delivers the 
exact asymptotics in the multi-qubit case (an earlier result is published in 
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(roland.hildebrandOimag.fr). This paper presents research results of the Action Concertée Incitative "Masses de 
données" of CNRS, France. The scientific responsibility rests with its author. 
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The paper is organized as follows. 

In the next section we characterize the extreme rays of the cone dual to the cone of ball-ball 
separable elements, namely the cone of linear maps that take the Lorentz cone to the Lorentz cone in 
the respective source and target spaces (Lorentz-to-Lorentz positive maps). It is well-known that the 
extreme rays of the dual cones characterize the largest faces in the primal cones JT]. These faces are 
of interest to us because they determine the radii of the largest separable balls around chosen elements 
in the separable cones. Note that describing the extreme rays of cones of Lorentz-to-Lorentz positive 
maps yields also a description of cones dual to ellipsoid-ellipsoid separable cones, because extreme rays 
are taken to extreme rays by invertible linear mappings. 

In Section 3 we describe the largest faces of ball-ball separable cones using the obtained families of 
extreme rays in the dual cones. This allows to get some insight into the structure of ball-ball separable 
cones. 

In Section 4 we compute the radius of the largest ball of ellipsoid-ellipsoid separable elements around 
the tensor product of points defining the central rays of the initial ellipsoid-generated cones. 

In Section 5 we apply this result to study the cone of separable unnormalized states of a multi-qubit 
system. We compute the above-mentioned lower bounds on the radii of largest separable balls around 
the unnormalized uniformly mixed state. 

In the last section we summarize our results. 



2 Extreme rays of cones of Lorentz-to-Lorentz positive maps 

In this section we compute the extreme rays of the cone of positive maps, i.e. those linear maps which 
take the Lorentz cone in the Lorentz cone in the respective source and target spaces. 
Let L m C R m , L n C R n be Standard Lorentz cones of dimensions m and n, i.e. 

L m = {(x ,xi,...,a; m _i) T |a;o > |(xi,...,x m _i) T |}, 
L n = {{ya,yi,...,y n -i) T \yo>\(yi,---,yn-i) T \}- 

We assume throughout the paper that min(n,m) > 2. Since L\ is isomorphic to the ray R+, the case 
min(n,m) = 1 is trivial. We call a linear map M : R m — * R" L m -to-L n positive or just positive if 
M[L m ] C L n . Since the Lorentz cones are self-dual, the cone V of such maps is dual to the cone of 
L m ® L„-separable elements. Moreover, as a consequence of this self-duality M is positive if and only 
if the adjoint map M T is positive in the sense that M T [L n ] C L m . 

In this section we determine the extreme rays of the cone V of positive maps. We represent maps 
from R m to R" bynxm matrices partitioned as 

where s is a scalar, h is a row vector, v is a column vector and A is a (n — 1) x (m — l)-matrix. Note 
that if M is a non-zero positive map, then the scalar s is strictly positive. 

Dcfinc two diagonal matrices J n = diag(l,I n -x), J m — diag(l, J m _i), where 1^ denotes the k x k 
idcntity matrix. Note that if x € R m is in the interior of L m , then x T J m x > 0. If x £ dL m , then 
x JfrjX - — ■ 0. 



Lemma 1 A map 



U =' t A ) eR " 



is L rn -to-L n positive if and only if 

3 X > : M T J n M h XJ m , \h\ < 1, 

or equivalently, 

3 X' > : MJ m M T t AJ„, \v\ < 1. 

Proof. By definition, M is positive if for any xo > 0, x G R" l_1 such that xq > \x\ we have yo > \y\, 
where 

' Vo \ _ M ( x o 



V 



We can rcwritc this equivalently as 



Vx ,x | x > 0, (x x T )J m ( ^° \ > : x a + hx > 0, (x x T )M T J n M ( X ° 



> 0. 
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By the 5-lemma 0,221 thïs is equivalent to the conditions 

(1 h) T £L m , 3 A > : M T J n M h XJ m , 

which gives the first set of conditions claimed by the lemma. The second set is obtained by considering 
the adjoint maps M T in the L n -to-L m positive cone. □ 

Let us establish necessary and sufficient conditions for a positive map to generate an extreme ray 
of the cone V . Note also that M generates an extreme ray if and only if M T generates an extreme ray. 
First we show that if a non-zero positive map M does not take the interior of L m in the interior of L n , 
then M is of rank 1. 

Lemma 2 Let M ^ be a positive map. Suppose there exists x 6 int L m such that Mx € dL n . Then 
the rank of M is equal to 1. 

Proof. Let the conditions of the lemma hold. Dcnote the point Mx £ dL n by y. Then y is 
contained in the linear subspace M [R m ] = ImM. Since M ^ 0, this subspace is non-zero. Moreover, 
since x £ intL m , there exists a neighbourhood U m C R m of x which is entirely contained in L m . Its 
image M[U m ] will be a neighbourhood U n of y relative to ImM. By the positivity of M the set U n 
is contained in L n and therefore in the face of y with respect to the cone L n . But y £ dL ni hence 
this face equals the intersection of L n with the linear subspace generated by y. Therefore y ^ and 
ImM = {ay \ a £ R}. Thus M has rank 1. □ 

Lemma 3 A positive map M of rank 1, partitioned as in £7J), generates an extreme ray if and only if 
\h\ = \v\=s. 

Proof. Let M be a rank 1 map, partitioned as in IjTJl. and let s > 0. Then we have 
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M=^ S v ) {sh) . 
The condition of positivity provided by Lemma ^ then takes the form 



(s h) T € L m , 3A>0: 1( ° T )(-<•')./„( j !>•/,).> \J„ 



s > \h\, 3X > : ^{s 2 - \v\ 2 ) ^ * T ^j(sh)h \J m 

^ s > \h\, s > \v\. (2) 

Hence a rank 1 map is positive if and only if s > \h\, s > \v\, s > 0. 

Let now M be a positive map of rank 1 and suppose that s > \h\ ^ 0. Then the maps 

M(A) = ~ ( * ) (s Xh) 

are positive for all A £ [—s/\h\,s/\h\] and M — M(l). The point A = 1 lics in the interior of the 
interval [— s/\h\, s/\h\], and the matrices M(A) are not múltiples of each other for different A. Hence 
M does not generate an extreme ray of V . 

A slightly modified argument can be applied if h = 0. Choose any vector h' with \h'\ = s and 
consider the maps 
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M(A) = - s I * j (s Xh') . 

Then M(A) is positive for all A £ [—1, 1] and M — M(0). Hence M cannot generate an extreme ray 
ncither. 

The same reasoning applies if s > \v\. Thus if M generates an extreme ray, then \h\ = \v\ = s. 

It rests to show that any rank 1 matrix with \h\ = |u| = s > generates an extreme ray. Let M 
be such a matrix. Suppose there exists a matrix SM such that M(X) — M + X5M £ V for all A in a 
neighbourhood U of zero. Let V be a neighbourhood of the unit vector eo £ R m that lies entirely in 
the interior of L m . Then for all A £ U and for all w £ V we have 

M(X)w = Mw + XSMw £ L n , Mw = -{s h)w ■ f * ^ £ dL n . 

Since A can vary in a neighbourhood of zero, the vectors Mw + X SM w have to lie in the face of Mw 
with respect to the cone L n . Then SMw lies in the tangent space to that face. This tangent space 
is the linear subspace generated by Mw. Hence SMw has to be a múltiple of the vector Mw. Since 
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this holds for all w £ V, the image of SM must be contained in the linear subspace of R™ generated 
by the set {Mw | w £ V}. Note that -(s h)w > 0, because w £ int L m . Therefore this subspace is 
one-dimensional and generated by (s v^) T . It follows that SM is of the form (s v T ) T u for some vector 
u £ R m . If we apply the samc linc of reasoning for the positive map M T , we conclude that 5M T is 
of the form (s h) T u' for some vector v! £ R n . Thus SM is proportional to M and M generates an 
extreme ray of V . 

This completes the proof of the lerrrma. □ 

It rcsts to consider the positive maps of rank strictly greater than 1. Let M be such a map, 
partitioned as in (JJJ. By Lemma|2|M takes the interior of L n to the interior of L m . Let Aut(L m ), 
Aut(L n ) be the automorphism groups of the cones L m , L n , respectively. We shall now show that if M 
generates an extreme ray, then there exist automorphisms U n £ Aut(L n ), U m £ Aut(L m ) such that 
U n MU m is doubly stochastic. (A positive map M is called doubly stochastic if M and M T take the 
central elements eo of the cones L m , L n into each other. Otherwise spoken, M is doubly stochastic if 
s = 1 and h = v = 0.) 

Define two functions p, q : R m — > R by p(x) = x T J m x, q(x) — x T M T J n Mx. Then the set 
N — {(p(x),q(x)) £ R 2 | x £ R m } is called the joint numerical range of the matrices J m , M T J n M 
underlying the quadratic forms p, q. It is known [2] that the set N is a convex cone. Lemma ^ states 
the existence of a number A > such that q(x) > Xp(x) for all x £ R m . Let A* be the maximal such 
A. Since M takes the interior of L m to the interior of L n , the set N has a non-empty intersection with 
the open first orthant. Therefore A* exists. Moreover, A* > 0, otherwise the matrix M T J n M would 
be positive semidefinite, which is not possible if the rank of M is strictly greater than 1. We have 
M T J n M - A* J m y and M T J n M - X* J rn - SJ m ^ for any d > 0. Hence there exists x ^ such 
that q(x) — \*p{x) = x T (M T J n M — X*J m )x = and p(x) = x T J m x > 0. Let us distinguish two cases. 

1. There exists a point x* £ R m such that q(x*) - X*p(x*) = and p(x*) > 0. 

Without restriction of generality we can choose x* such that x* £ intL m and p(x*) = 1. Denote 
Mx* by y*. Since M takes intL m to intL n , we have y* £ intL n . In fact, (y*) T J n y* = q(x*) = X* > 0. 
Let U m , U n be automorphisms of the cones L m , L n , respectively, preserving the quadratic forms J m , 
J„, respectively, such that U m x* — e™ £ R m and U n y* = V~X*eQ £ R". (Here e™,eg are the unit 
vectors in the direction of the coordinate xq in the respective spaces.) Such automorphisms exist since 
the Lorentz cones are homogeneous |15| . 

Define a map M = (X*)^ 1 ^ 2 U n MU rn 1 . By the positivity of M this map is also positive. We have 
M e ™ = {X*Y 1 l 2 V n MV^U m x* = (A*)- 1 / 2 ^^* = e#. Since x* is contained in the nullspace of the 
positive semidefinite matrix M T J n M — X*J m , we have M T J n Mx* = M T J n y* = X* J m x* . It follows 
that 



M T e% = M T (J n e%) = {X*)-^U~ T M T U^{U- T J n U?)<% - {X*)- 1/2 U m T M T J„[(A*)" ^y*} 
= (A*) 1 U m T [X* J m x*] — U m T J m x* — J m U m x* — e™. 



Hence M is doubly stochastic. 

Remark: A similar statement for cones of maps that take the positive semidefinite cone to the 
positive semidefinite cone was proven by Leonid Gurvits |3] . 

2. For any point x £ intL m we have q(x) > X*p(x). 

We noted above that there exists x* ^ such that q(x*) — X*p(x*) and p(x*) > 0. Since p(x) > 
yields q(x) > X*p(x), we have p(x*) — 0. Without restriction of generality we can choose x* such 
that x* £ dL m . Denote by Ln the nullspace of the positive semidefinite matrix M J n M — X* J m . 
This linear subspace contains x* £ dL m and does not intersect the interior of L m . Hence it lies in the 
orthogonal complement to the element J m x* £ dL m . 

On the other hand, (M T J„M - X*J m )x* = yields for any vector v £ KerM that v T (M T J n M - 
X*J m )x* — — X*v T J m x* = 0. Hence the kernel of M lies also in the orthogonal complement of J m x* ■ 
Equivalently, J m x* lies in the image of the matrix M T and there exists a vector v £ R™ such that 
J m x* = M T v. 

Let now A = J n v(J m x*) T and consider the family of maps M(a) = M + aA. For any vector 
w £ Ln and for any a we have 



because Aw = J n v(J m x* ,w) = and A T J n Mw = J m x*v T J% J n Mw = J m x*(J m x*,w) = 0. Hence 
there exists 5 > such that for all a £ (—6, +S) the matrix M(a) T J n M(a) — X*J m lies in the face of 
the positive semidefinite cone generated by the matrix M T J n M — X*J m . 
Let the matrix M(a) be partitioned as 



[M{a) T J n M{a) ~ X* J m ]w 



[aA T J n M + aM T J n A + a 2 A T J n A]w = 0, 
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We have (s(0) h(0)) — (s h) £ intL m , because otherwise the positive map M T would take the vector 
e 7 Q £ ínt L n to the vector (s h) T £ dL m , and M T would have rank 1 by Lemma [21 Therefore there 
cxists 6' > such that (s(a) h(a)) T £ int L m for all a £ (—6', +6'). 

Then by Lemmanthe map M(a) is positive for all a with |a| < min(#, 5') and hence contained in 
the cone V. Since the rank of A equals 1, but the rank of M is strictly greater than 1, the matrices 
M, A cannot be collinear. It follows that M does not generate an extreme ray of V. 

We have proven the following 

Corollary 1 Let M be a positive map of rank strictly greater than 1 and let M generate an extreme 
ray ofV. Then there exist automorphisms U n £ Aut{L n ), U m £ Aut(L m ) such that U n MU m is doubly 
stochastic. □ 

Note that for any automorphisms U n £ Aut(L n ), U m £ Aut(L m ) the matrix U n MU m generates an 
extreme ray of V if and only if M generates an extreme ray of V . Let us characterize the extreme rays 
that are generated by doubly stochastic matrices. From Lemma^it follows that a doubly stochastic 
matrix, partitioned as in l|T|l. is positive if and only if cr max (^4) < 1, where <r max denotes the maximal 
singular value. 

Lemma 4 Let M be a doubly stochastic positive map, partitioned as in JT\). and let M generate an 
extreme ray ofV. Then all singular vàlues of A equal 1. 

Proof. Let us assume the contrary. Suppose M is doubly stochastic and positive, partitioned 
as in QJ, with a m i n (A) = ò < 1. Let A = UDV be the singular value decomposition of A and 
oi, . . . , c m in(m-i,n— i) its singular vàlues in decreasing order. Here U, V are orthogonal matrices of 
appropriate size and D = diag{a~x, 02, . . . , 0" m in(m— i,n-i)) is a (n — 1) x (m — 1) matrix with the 
singular vàlues of A on its main diagonal, all other elements being zero. Note that 0'min(m-i,n— 1) = 
<7 < 1. Let us dcfinc an afhne one-parametric family of diagonal (n — 1) X (m — 1) matrices by 
D(a) — diag{a\, 02, . . . , 0'min(n,m)-2) a )- Then the maps 



M(a) 



1 

UD{a)V 



are positive and hence belong to V for all a £ [— 1, +1]. Note that M — M(a). Since these matrices are 
not proportional for different vàlues of a, and £ (— 1, 1), the map M does not generate an extreme 
ray of the cone V. This proves the lemma. □ 

Lemma 5 Let M be a doubly stochastic positive map, partitioned as in yj), and let all singular vàlues 
of A equal 1. Then M generates an extreme ray ofV if and only í/min(n, m) > 2. 

Proof. Let M be a map satisfying the assumptions of the lemma. Assume also without restriction 
of generality that n > m. Then we have A T A — I m -i- Let us first show that M does not generate an 
extreme ray if min(n, m) = 2. 

If m = 2, then the matrix A is a unit length column vector. Consider the two maps 

Mi = ( \ \ (1 1), M 2 = ( } A ) (1 - 1). 

These maps are positive by condition J2Jl and not proportional. Moreover, we have M = \ {M\ + M2). 
Hence M does not generate an extreme ray of V . 

Suppose now that n > m > 3. Assume there exists an n X m matrix 



h s 
vs A s 



and a number e > such that the map M(a) = Al + aM$ is positive for all a £ (— e, +e). The 
assumption that the upper left element of M$ is zero does not restrict the generality, because this 
element can be made zero by adding to Mg an appropriate múltiple of M. Let us develop the positivity 
condition of Lcmman] We have that M(a) is positive if and only if \a\\hs \ < 1 and there exists A > 
such that 



1 - A - a 2 \vs\ 2 a(h 5 -vjA) - a 2 vjA s 

a(hj - A T v 5 ) - a 2 Ajv 5 (A - l)J ro _i - a(AjA + A T A S ) + a 2 (hjh 5 - AjA 5 ) 



h 0.(3) 



Here A may depend on a. We obtain in particular —a(AjA+A T As)+a 2 (hJhs—AjAs) >z (1— A)J m _i >z 
>Z 0. A necessary condition for this inequality to hold for all a £ (— e, +e) is that 
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AjA + A T A S = 0. It follows that hjh s h AjA s + \v s \ 2 I m -i- The left-hand side of this inequality 
is a matrix of rank not exceeding 1, while the right-hand side is positive semidefinite. Hence the 
rank of the right-hand side cannot exceed 1 too. Since m — 1 > 2, it follows that vg — and 
As is of the form whs, where w is a column vector of appropriate size. This yields the inequality 
a 2 hjhs(l — \w\ 2 ) >z (1 — A)7 m _i, which implies A = 1 for a similar reason. But then the upper left 
element of matrix is zero. Therefore a(hg — vjA) — a 2 vjAs = ahs = for all a G (— e, +e) and 
hs = 0, As = whs = 0. This proves that M generates an extreme ray of V. □ 

Combining the results obtained so far, we can characterize the extreme rays of the cone V as follows. 

Lemma 6 Let the positive map M be partitioned as in £7J) and suppose that it generates an extreme 
ray of the cone V . Then either M is of rank 1, with \h\ — \v\ — s, or there exist automorphisms 
U m G Aut(L m ), U n G Aut(L n ) such that 

U n MU m = ( J °, 

with all singular vàlues of A' equal to 1 . 

If min(r7i, n) > 3, then all matrices of the above types generate extreme rays. If min(r7i, n) = 2, 
then only those of them which are of rank 1 generate extreme rays. □ 

Note that for any pair of non-zero elements x, y G dL m in the boundary of the Lorentz cone L m 
there exists an automorphism U m of L m that takes x to y. Further, for any orthogonal matrix C/ m _i 
of dimension (m — 1) x (m — 1) the m x m matrix diag(l, Í7 m _i) represents an automorphism of L m . 
This allows us to reduce the extreme rays of V to two canonical forms. Define the two positive maps 

M^diagíí ! ! Vo,...,oW n l2x2 n ° 2 *^ 



11/ / \ 0(n-2)x2 0(n-2)x(m-2) 

M 2 = diag(l, . . . , 1) = <( V. (n _ m)xm )' n ~ m ' (4) 

(In 0„ x(m _„)), n < m. 

Here lfcxZjOfcxí denote k x l matrices filled with ones and zeros, respectively. 

Definition 1 We call a positive map M of Type I if there exist automorphisms U m G Aut(L m ), U n G 
Aut{L n ) such that U n MU m = M\. We call M of Type II if there exist automorphisms U m G Aut(L m ), 
U n G Aut(L n ) such that U„MU m = M 2 . 
We have the following theorem. 

Theorem 1 Let the positive map M generate an extreme ray of the cone V . Then M is either of Type 
I or of Type II. 

If min(m,rí) > 3 7 then all matrices of Types I and II generate extreme rays. 7/min(rn,n) = 2, then 
only the matrices of Type I generate extreme rays. □ 

The theorem shows that the structure of the cone V of positive maps is more complex than the 
structure of the Lorentz cones Lk, but is still relatively simple. While the Lorentz cone has only one 
kind of extreme rays (which are equivalent with respect to the action of the automorphism group), 
the cone of positive maps has two kinds. An exception are the cones of L^-to-Ln positive maps. In 
this case the extreme rays form two copies of the boundary dL n of the cone L n which are located in 
mutually orthogonal subspaces. 



3 Largest faces of ball-ball separable cones 

In this section we give a description of the largest faces of ball-ball separable cones, departing from the 
two families of extreme rays of the cone of positive maps obtained in the previous section. 

We call an element B of the space R m <g)R n L m ®L n -separable or just separable if B can be expressed 
as a finite sum J2k=i Xk ® V k °f P r °duct elements such that Xk G L m , yk G L n for all k = 1, . . . , N. The 
separable elements form a convex cone in R m <g>R™, the separable cone, which will be denoted by K sep . 
This cone is dual to the cone V of positive maps considered in the previous section. For convenience 
we will represent the elements of R m % R™ asnxm matrices such that the scalar product of a linear 
map M : R m — > R" with an element B G R™ ® R" is given by (B, M) = tr (M T B) = tr (B T M). In 
this representation a product element x <g> y is given by the rank 1 matrix yx T . 

It is well-known that the largest faces (i.e. non-trivial faces that are not an intersection of other, 
strictly larger faces) of a convex cone K have the form {x G K \ (x,y) — 0}, where y generates an 
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extreme ray of the dual cone K* jllj . Let us compute the faces corresponding to the extreme rays of 
the cone of positive maps V described by Tlieorem^ By this theorem, there are two kinds of extreme 
rays. These generate two kinds of largest faces of the separable cone. 
Let us define two Standard faces of the separable cone K sep by 

Fi = {B £ K sep | (B, Mi) = 0}, F 2 = {B £ K sep \ (B, M 2 ) = 0}, 

where M\ , M 2 are the positive maps (J2J . 

Definition 2 We call a face F of K sep of Type I if there exist automorphisms U m £ Aut(L m ), 
U n £ Aut(L n ) such that {U n BU m | B £ F} = F\. We call a face F of K sep of Type II if there exist 
automorphisms U m £ Aut(L m ), U n £ Aut(L n ) such that {U n BU m \ B £ F} = F 2 (or U n FU m = F 2 
for short). 

Hence all faces of Type I are affinely isomorphic to Fi, while all faces of Type II are affinely 
isomorphic to F 2 . Let us determine the structure of these two sets. 

Proposition 1 Fi is affinely isomorphic to the convex conic hull of the union 



Zi 



z — ( z 0j %ll ■ • ■ ; Z n+m ^ 2 ) ri 



1,(Z!-1) 2 +^Z2 = 15 



z. 



n+m— 2 



o \ u 



fc=2 



U < Z = (z , Zi, . . . , z„+ m -2) J 



n+m— 2 



z = 1, Zi 



_x = 0, (z m - l) 2 + *k = 1 \ c R 



n+m— 1 



k—m+1 



Remark: Thus a section of the cone F\ is affinely isomorphic to the convex hull of two intersecting 
spheres S™ -1 , S n ~ 1 which are located in orthogonal subspaces. 

Proof. The set F\ = {B £ K sep \ (B, M\) = 0} is given by the convex hull of those extreme rays 
of K sep that are orthogonal to M\. The extreme rays of K sep are tensor products of the extreme rays 
generating the individual Lorentz cones L m , L ni i.e. generated by elements of the form 



B = 



1 

h T 

> n—X 



1 h 

v vh 



where h £ R m is a row vector, v £ R™ is a column vector with \h\ — \v\ — 1. We have 



(5) 



(B,M 1 )=tr 



1 h 

v vh 



(eT'Y 



Jfi— 1 n—1 ( m— 1 \T 



5 l - 1 (e™- i ) J 



= (l + (v,e n r 1 ))(l + (h,e?- 1 )). 



Here e\ is the unit vector in the direction of the first coordinate in the space R fc . Therefore (B, M\) = 

n-l 

" H 



if and only if v = — e™ 1 or h 



-eT 1 . Thus we obtain 



Fi = conv 



1 



x J ' + y(l -(ef-T) 



x £ dL m , y £ dL n 



(6) 



= < 



' / x + y xt-yo x 2 
-xq + yx —xi - yi -x 2 
V2 -yi 



*Em— 1 ' 
*^m— 1 



o 



(n— 2)x(m— 2) 



k \ 2/n-l -2/n-l 

It is now easily seen that the affine map / : R, n + m - 1 



a; > |(x l5 . 

yo > \(yi, 



) 



R" 



R™ given by 



i %m— l) Is 
^n-l) T | 



> . 



2f) 



•^n+m — 2 



( z zi-l 

Z m — 1 2 — Zo — Zl — Zn 
Zm+1 Z m ^i 



Z2 
-2 2 







(n— 2)x(m— 2) 



\ ^n+m-2 



■^n+m— 2 



Zm-1 \ 
Zm—1 



) 



is an affine bijection between the union Z\ and a set of generators of the cone Fi. □ 

Let us now consider the second kind of largest faces. Denote by Sym(k) the space of real symmetric 
k x k matrices. 
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Proposition 2 The face F 2 is affinely isomorphic to the set 

Z 2 = {Ae Sym(mía{n, m)) \ A > 0, A Q0 = tr A/2}, 

where Aqq is the upper left element of the matrix A. 

Proof. Let n > m without restriction of generality. Then the positive map 



1 Olx(m-l) 
M 2 = | 0( m _!)xl —Im-X | = M 2 

0(n— m)xl 0(n— m) X (m— 1) 



1 

'(m-l)xl 







lx(m-l) 
-Im-1 



generates an extreme ray of V and is of Type II. Instead of the face F 2 we will consider the isomorphic 
face F 2 = {B G K sep \ (B, M 2 ) = 0}. 

This face is given by the convex hull of those extreme rays of K sep that are orthogonal to M 2 . 
Let such an extreme ray be generated by the tensor product 10 . Let the vector v be partitioned in a 
subvector v a of dimension m — 1 and a subvector vi, of dimcnsion n — m. We have 



(B, M 2 ) = tr 



1 h 

v vh 



Mo 



= tr 



1 

h T 



-h T v? 



1 - hv a - 



Note that \h\ — 1, |« | < 1. Therefore (B, M2) = if and only if v a = h T and = 0. Thus F2 is given 
by the convex conic hull of the set 



1 

h T 



h 

h T h 



0(n— m)xl "(n m)x(m 1) 



This hull is equal to the set 



.4 







(n— m) xm 



o 



.4 

(n— m) x m 



A e Sym(m), AyQ,rkA<l, A oa = trA/2 



A e Sym(m), A^O, A oa = tr A/2 



The last relation is a consequence of the following fact. 

// L is an linear subspace of Sym(k) of codimension 1, then the intersection of L with the cone 
S+(k) of PSD matrices in Sym(k) equals the convex conic hull of all rank 1 PSD matrices which are 
contained in L. 

Indeed, if these two sets do not coincide, then there exists a linear functional L 1 on Sym(k) that 
strictly separates some point in In S+(k) from the convex conic hull of all rank 1 PSD matrices in 
L. But this contradicts the convexity of the joint numerical range 2 of the quadratic forms on R fc 
induced by L and l·I . 

This completes the proof. □ 

Above characterizations of the Standard faces F\ , F 2 allow us to characterize all faces of Types I 
and II. 

Lemma 7 The faces of Type I are parameterized by a pair of vectors (h, v), where h G S"™~ 2 C R' Tl_1 
is a row vector of length 1 and v S S n ~ 2 C R" _1 is a column vector of length 1. The face Fj(h,v) 
corresponding to such a pair (h, v) is given by 



Fz(h, v) = conv 



1 



2/(1 h) 



x e dL m , y G dL n 



Proof. Let F be a face of Type I. Then there exist automorphisms U n , U m of L n , L m , respectively, 
such that F = U n F\XJ m . Define the vectors 

tí = Un ( ) e dL n , v' = Ul(^ _ e Li ) G 3L rn . 

Let h = (h') T /h' Q , v — v'/v' be normalized múltiples of h',v'. Then description of the Standard 
face Fi shows that F has the form Fj(h,v) defined in the theorem. 

On the other hand, the generator sets of Fj(h, v), Fx(h', v') are different whenever (h, v) ^ (h 1 , v'). 
Hence F T (h,v) ^ F^h' ,v') for (h,v) ± (h',v'). □ 



s 



Corollary 2 Any two faces of Type I have a non-trivial intersection. 

Proof. Let Fi(h, v), Fj(h', v') be two faces of Type I. Then the elements 

( v' v'h^)'(v vtí ) G Ksep 
are contained in both Fj(h,v) and Fi(h',v'). □ 

Corollary 3 Any face of Type I has a non-trivial intersection with any face of Type II. 

Proof. Let Fi(h,v) be a face of Type I, and let Fjj be a face of Type II. Then there exist auto- 
morphisms U n ,U m such that Fjj = U n F2U m . We assume n > m without loss of generality. Choose 
v' £ dL n such that U- l {l (v') T ) T {l WjU' 1 is in F 2 . Then the element (1 (v') T ) T (í h) is shared by 
the faces F H and Fj(h,v). □ 

On the other hand, faces of Type II do not necessarily have a non-trivial intersection. 

In this section we have described the largest faces of the separable cone K sep . There are two types 
of such faces. All faces of one type are equivalent with respect to the action of the automorphism 
groups of the underlying Lorentz cones. Any other non-trivial face is an intersection of some largest 
faces. The faces of Type I are affinely isomorphic to the convex conic hull of two spheres S n ~ 2 ,S m ~ 2 
which intcrsect each other in onc point, but lic in orthogonal subspaces. The faces of Type II are 
intcrsections of the cone of positive semidefinite min(n, m) x min(ra, m)-matrices with a linear subspace 
of codimension 1. Note that the manifold formed by the union of relative interiors of Type I faces 
has (n + m — 1) + (n — 2) + (m — 2) = 2(n + m) — 5 dimensions, whereas the boundary of K sep is 
nm — 1-dinicnsional. But 2(n + m) — 5 < nm — 1 if min(n, m) > 3. Hcnce the boundary of the 
L2 <B> L„-separablc cones is formed by Type I faces, while the boundary of K sep for min(n, m) > 3 is 
formed by Type II faces. 



4 Radii of largest separable balls 

Extreme rays and largest faces remain invariant under linear bijections. Therefore the results obtaincd 
in the last two sections are extendible to cones that are separable with respect to linear images of 
Standard Lorentz cones. In this section we compute radii of largest separable balls. These radii are 
naturally invariant only under orthogonal mappings, therefore results obtained for L m ® L m -separable 
cones will not extend to arbitrary linear images of the Lorentz cones. In order to cover this more general 
case, we will consider general ellipsoid-generated cones. We shall compute the radius of the maximal 
ellipsoid-ellipsoid separable ball around the tensor product of elements generating the central rays of the 
two individual ellipsoid-generated cones. Let us first give a precise definition of an ellipsoid-generated 
cone and its central ray. 

Let B C E be a closed sòlid ellipsoid with nonempty interior in some n-dimensional real vector 
space. Suppose that the origin of the space is not contained in B. Then the conic hull of B is the 
image of the Standard Lorentz cone L n under a regular linear mapping. Moreover, by a rotation it can 
be transformed to some standardized cllipsoidal conc 

Kst(P) = j(x ,2!i, ■ • ■ ,a;™-i) T | xo > Vx T Px, x = (xi, . . . ,x„_i) T | , 

where P is a positive definite symmctric [n— 1) x (n— l)-matrix. The set of positive defmite symmetric 
(n — 1) x (n — l)-matrices parameterizes the set of standardized ellipsoidal cones in R". Wc define the 
central ray of K st {P) as the ray generated by the unit vector eo = (1, 0, . . . , 0) T . 

Let now K\ C R m , K2 C R" be standardized ellipsoidal cones given by positive definite matrices 
Pi, P2 of dimensions (m — 1) x (m — 1) and (n — 1) x (n — 1), respectively: 

K 1 = K st (Pi) = |(x ,xi, . . . ,a; m _i) T € R m | xq > y / x T P 1 x 1 x = (x u . . . ,a; m _i) T | , 

K 2 = K st (P 2 ) = {(y ,yi,...,y n -i) T € R" | yo > \/y T P2V, y = (yi, ■ ■ ■ ,2M-i) T } • 

Dcnote by e^, e^, . . . , e^ f _ 1 the unit vectors along the coordinate axes of the space R N . Then the unit 
vectors e™, define the central rays of the cones Ki,K 2 . 

As in the previous section, we shall call an element B £ R m (g)R rl Ki(g)K 2 -separable or just separable 
if B can be expressed as a finite sum J2k=i Xk ® V k °^ P r °duct elements such that x^ £ K\ , y^ £ K 2 
for all k — 1, . . . , N. The cone of K\ ® inseparable elements, the separable cone, will be denoted 
by K sep . We will represent the elements of R m ® R" as n x m matrices such that a product element 
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x ® y is given by the rank 1 matrix yx T . Then the product e™ 
elements everywhere except a 1 in the upper left còrner. 



i is given by a matrix that has zero 



We shall compute the radius of the largest ball around the unit length vector e™ 



£ R" 



R" 



consisting of K\ ® inseparable elements. A ball B C R m ® R" 
only if the conical hull of B is contained in the separable cone K sep 
cone, and as such an ellipsoid-generated cone. 



consists of separable elements if and 
. This conical hull is a ball-generated 



Lemma 8 Consider the real vector space R . The cone generated by a ball of radius p < 1 around 
the unit vector eo equals the standardized ellipsoidal cone K s t{r~ 2 1 n -i) with r = . p 

Proof. By the rotational symmetry of the ball-generated cone it must equal a standardized ellip- 
soidal cone K st (P) with the matrix P being proportional to the identity matrix. From the definition 
of K s t{r~ 2 1 at -i) it follows that r is the radius of the ball created by the intersection of the cone 
Kst(r~ 2 In-i) with the hyperplane given by the equation xq — 1. The relation r 



consequence of the similarity of appropriate rectangular triangles formed in the eo-ei plane of R w . □ 
Note that r = p is a monotonous function of p £ [0, 1). 

Let us denote the ball-generated cone K s t(r~ 2 I nm -i) by Kb a ii(r)- Identify R Iim with R m <g)R™ by 
idcntifying the basis vectors e^n+V k = 0, . . . ,m — 1, l = 0, . . . ,n — 1, with the orthonormal basis of 
tensor products e™ <X> ef. Then the cone Kb a ii(r) is generated by a ball centered on e™ <g> efi. 

Let K*,Kbail(r)* denote the cones dual to K sep , Kb a ii(r). These cones reside in the space 
(R nm )* = R m ™, whose elements will likewise be represented by n x to matrices. The scalar product of 
a matrix B £ R m ®R™ with a matrix M £ (R" m )* will be defined as (B, M) = tr (M T B) = tr (B T M). 

Lemma 9 Let r be the largest number such that the inclusion K^aii (r) C K sep holds. Then 



max \ hP x h T + v T P 2 v + tr (A T P 2 APi) 



M 



h 
A 



is L rn -to-L n positive 



(7) 



Proof. We have Kb a n(r) C K sep if and only if K* ep C Kball ( r )*- By means of Standard linear 
àlgebra one establishes that K ba u{ r )* — ^a«( r_1 ), K st {P)* — K st (P^ 1 ) and K* ep is the cone of 
Ki-to-K^, or K st (Pi)-to-K st(-P 2 _1 ) positive maps. It follows that the largest r satisfying the inclusion 
Kbaii{r) C K sep equals the inverse of the smallest R such that any ií sí (Pi)-to-fí st (P 2 " 1 ) positive map 
lies in the cone Kb a u(R)- 

Let us characterize the cone of K s t(Pi)-to-K st (P2 1 ) positive maps and the cone Kb a ii(R). Let 
M : R m — * R" be a linear map, partitioned as 



M = 



1 h 

v A 



(8) 



where h is a row vector of length to — 1, v is a column vector of length n— 1 and A is a (n — 1) x (m — 1) 
matrix. 

Since for a positive definite matrix P and for any vector x we have V x T Px = \P x / 2 x\, we can 
characterize the cones K%, K~2 as follows: 



Kl = 



|(x ,a;i, . . . ,x m _i) T | (x , (xx, . . . , x m -x)Pl /2 ) T £ L m \ , 

{(2/0,2/1, ■ ■ ■ ,2M-l) T | (2/0, (2/1, ■ ■ ·,J/n-l)-P2~ 1/2 ) T e L n\ ■ 



It follows that the map diag(l, P 1 1 ^ 2 ) is an isomorphism between L rn and K\ and the map diag(l, P 2 1 ^ 2 ) 
an isomorphism between if| and L n . Hence M is Kx-to-K^ positive if and only if the map 



M = diag(l, P 2 1/2 ) • M ■ diag(l } P 1 = 



-l/2s _ 



.4 



Pn 



-1/2 



V P n 



hP- 

-1/2 



-1/2 
L 

AR 



-1/2 



is L m -to-L n positive. 

Let us examine the cone Kbaii(R)- By definition M is in Kbaii(R) if Mqo > ^Jj2(k 1)^(0 0) M^/R 
(here M^; denotes the elements of M). If M is partitioned as in (JSJ, then M o = 1 and M G Kb a ii{R) 



if and only if R > . /j^ 



/(fc,0^(o,o) ^fei 



ll^ll 
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Hence we obtain the following characterization of the largest number r such that any K\-to-K 2 
positive map is contained in K^ai^r^ 1 ): 

1 hp- 1/2 ^ ^ 



r- 1 = max{ J\h\ 2 + \v\* + \\A\\l 



max a/I^i 72 ! 2 + l^l 2 + WP^AP^Wl 



M=\ p -i/2 v p -i/2 Ap -i/2 J is L m -to-L n positive j 

^ ^ is L m -to-L n positive j- 



max <^ hP x h T + v T P 2 v + tr {A T P 2 AP 1 ) 



M = ( ^ ^ j i> i,,,,-: <»-!,„ [)i»iti\c y L 



We shall now calculate expression Q. We have to compute the maximum of the function F(M) = 
hPih T + v T P 2 v + tr (A T P 2 APi) over the set S of L m -to-L n positive maps M which are partitioned as 
in (JHJ, i.e. with the upper left element being equal to 1. We shall show that F achieves its maximum 
either at a rank 1 map or at a doubly stochastic map. 

The following lemma is verified by direct calculation. 

Lemma 10 For any integer n > 2 and any row vector b € R™ -1 the linear transformation 

I 1+|fc|2 H J "-ï-(TTW; 

U n (b) 




b T b \ 1 
L (Ï+Ï6FF ) 

preserves the quadratic form J n , i.e. J n = U n (b) T J n U n (b) , and is hence an automorphism of the cone 
L n . □ 

Let M be an L m -to-L n positive map, partitioned as in JSJ. By the preceding lemma the maps 
U n (b n )M, MU m (b m ) are also positive for all b n S R n_1 , b m S R m_1 . The upper left elements of these 
products are given by 

rni(b n ) = 1 + \b n \ 2 + - \ n 1 ^ ^p^- v > 0, m r (b m ) = 1 + \b m \ 2 + hbj n > 0. 

Consider the famílies of positive maps M;(6 n ) = U n (b n )M/mi(b n ), M r (b m ) — MU m {b m )/m r (b m ), 
parameterized by row vectors b n E R" _1 , b m G R m_1 . The upper left element of the corresponding 
matrices equals 1. Hence M/(6 n ), M r (b m ) can be partitioned as in JSJ: 

" ( MK) A&) ) • - ( mL) A lt) ) • <« 

where hi(b n ), vi(b n ), Ai(b n ), h r {b m ), v r (b m ), A r {b m ) are vectors and matrices dcpcnding accordingly 
on the parameter vectors b ni b m . Define two scalar functions 

Fi{b n ) = F(Mi(b n )) = hi(bn)PLhi(K) T + vi{b n ) T P 2 vi{b n ) + tr (^(6„) T P 2 A i (6„)P 1 ), 
F r (b m ) = F(M r (b m )) = hrib^Pihribmf + v r (b m ) T P 2 v r (b m ) + tr {A r (b m ) T P 2 A r (b m )P 1 ). 

Lemma 11 Let a L m -to-L n positive map M , partitioned as in |3J), realize the maximum of the function 
F. Then M generates an extreme ray of the cone V of L m -to-L n positive maps. The corresponding 
functions Fi(b n ), F r (b m ) have global màxima at b n = 0, b m = 0, respectively. As a consequence, their 
gradients at b n = and b m = vanish. 

Proof. The function F(M) = hP\h T + v T P 2 v + tr (A T P 2 APi) is strictly convex on the convex set 
S. Hence its maximum is achieved at an extreme point of this set. Equivalently, the map realizing the 
maximum of F generates an extreme ray of V . 

Let the map M realize the maximum of F. Define the families of maps JÜJ). Now note that 
U n (0),U m (0) are the identity maps, hence M;(0) = M r (0) — M. Since the maps M/(ò„), M r (b m ) are 
in S for all b n ,b m , the functions Fi(b n ), F r (b m ) attain their global màxima at the origin. □ 

Lemma 12 Let an L m -to-L n positive map M , partitioned as in f%[). realize the maximum of F. If 
M has rank 1, then this maximum is given by F max = —1 + (1 + A max (í'i))(l + A max (P2)) (here A ma x 
denotes the maximal eigenvalue). 
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Proof. Let Al satisfy the assumptions of the lemma. If M is of rank 1, then A = vh. By the 
previous lemma M generates an extreme ray of V . By Lemma |3]we then have \h\ = \v\ = 1. On the 
other hand, any pair of unit length vectors {hi , v') defines a positive map of rank 1 via 



We have 



F(M(h',v')) = h'PrQi'f + {v') t Pív' + {v') T P 2 v' + {v') T Piv'tíP l {h') T 

= -\ + {l + h l P 1 {h l ) T ){\ + {v') T P 2 v l ). 

It follows that 

F max = max F(M(ti,v')) = -1 + (1 + A max (Pi))(l + A max (P 2 )). □ 

Lemma 13 Let M be an L m -to-L n positive map, partitioned as in @). Then <J max (A) < 1. 

Proof. Let M be a map satisfying the assumptions of the lemma and suppose that a — er max (A) > 1 . 
Then there exist unit length column vectors u, w of appropriate dimensions such that au — Aw and 
hence a — u T Aw. Without restriction of generality we can assume that hw — u T v < (otherwise we 
multiply u, w by —1). Then we have 



i y i h \ i 

— u I \ v A I \ w 



But 

1 \ fi 



= 1 + hw — u v — er<l — er<0. 



, e L n , e L m , M I e L n , 

— u / \ w J \ w J 

the last inclusion being due to the positivity of M . Hence the scalar product of two vectors in L n is 
negative, which leads to a contradiction with the self-duality of L n . Thus the assumption cJmax(^4) > 1 
was false, which completes the proof. □ 

Lemma 14 Let an L m -to-L n positive map M, partitioned as in realize the maximum of F. Sup- 
pose further that this maximum is strictly greater than the maximum over the rank 1 maps established 
in Lemma, Then h = v = 0. 

Proof. Let M satisfy the assumptions of the lemma. We shall now compute the gradients of the 
functions Fi(b n ), F r (b m ) at b n = 0, b m = 0. We have 

U' n (b) | b =o = y ( ò Ty o)' m 'i( b )\b=o = b'v, m' r (b)\ b = = h(b T )' . 

Hence we obtain 

fcí(Mk=o = b> n (A - vh), «í(Mk=o = (In-i - vv T ){bl)', Mb n )\ bti=0 = (bl)'h - {v T {bl)')A, 
K(b m )\ bm=0 = b' m {I m - X - h T h), <(6 ro )| 6m=0 = (A- vh){b T m )', A r (b m )\ bm=0 = vb' m - A(b' m h T ). 
It follows that 

P/(ò«)k=o - 2(^Pi/i T + {vJ)'P 2 v + tr{{Aj)'P 2 AP l )) 

= 2b' n [(/„_! + P 2 )AP x h T + {-{hPxh T + v T P 2 v + tr(A T P^P^In^ + P 2 )v] , 
F;(6 m )| hm =o = l(h! r P x h T + (vJ)'P 2 v + tr((Aj)'P 2 AP 1 )) 

= 2b' m [(J m _i + P 1 )A T P 2 v + (-(hPih T + v T P 2 v + tr(A T P 2 AP 1 ))I m - 1 + Px)h T ] . 

Sincc the vanishing of the gradient is a necessary condition of maximality of the functions Fi,F r , we 
obtain the equations 

(/„_! + P 2 )AP x h T = [(1 + hP x h T + v T P 2 v + ír(A T P 2 ylP 1 ))/„_ 1 - (/„_! + P 2 )]v, 
(7 m _i + Pi)A T P 2 w - [(1 + hP x h T + v T P 2 v + ír(A T P 2 ylP 1 ))/ m _ 1 - (7 m _! + P)]/i T . 

The maximum of P is given by P max = 1 + hP x h T + v T P 2 v + tr(A T P 2 AP 1 ). It follows that 

AP x h T = [(J n _i + P^^P^ - / n _iK A T P 2 v = [(7 m _! + Pi)"^^ - 7 m _i]/i T . (10) 
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By the assumptions of the lemma we have P max > (1 + A max (Pi))(l + A max (P 2 )). Therefore {I n -i + 
P2)~ 1 F max - I n -x y A max (Pi)J„_i and (J m _i + Pi) _1 P m ax - I m -i >■ A max (P 2 )/ m _i. Hence the 
matrices on the right-hand sides of I|1U|) are invertible and h = implics v = and vice versa. Let us 
assume that h ^ and u ^ 0. Taking the norms on both sides of equations (JTÜJ, we get 

|H| 00 A n , Mt (Pl)|/l| > |APl/l T | = |[(/ n _i + P 2 )" 1 Pmax " /„-l]«| > A m „(i\)|«|, 

P|| co A max (P 2 )|i;| > \A T P 2 v\ = |[(/ m _! +Pi)- 1 P max -/ m -i]/i T | > A max (P 2 )H. 

Combining, we obtain ||-<4.||oo = ffmax(^) > lj which by Lemma |13I leads to a contradiction with the 
positivity of M. Hence h = v = 0, which completes the proof. □ 

The lemma implies that if M realizes the maximum of F and has a rank greater than 1, then it 
must be doubly stochastic. 

Lemma 15 Let Ai(Pi), A 2 (Pi), . . . , A m _i(Pi) and Ai(P 2 ), A 2 (P 2 ), . . . , A„_i(P 2 ) be the eigenvalues of 
the matrices Pi, P 2; respectively, in decreasing order. Then the maximum of F is given by the expression 

P max = max{-l + (1 + Ai(Pi))(l + Ai(P 2 )), Er=i ( "' mKl A fc (Pi)A fc (P 2 )}. 

Proof. We have shown above that the maximum of F is achieved either at a rank 1 map, in which 
case it equals F max = — 1 + (1 + Ai(Pi))(l + Ai(P 2 )), or at a doubly stochastic map. Suppose we are 
in the second case, and the map rcalizing the maximum of F is partitioned as in JSJ with h = v = 0. 

Since the maximum is achieved at a map generating an extreme ray of the cone V, all singular 
vàlues of the matrix A equal 1 by Lemma 0] Assume without restriction of generality that n > m. 
Then the singular value decomposition of A is given by 

A = UDV = u( Im - 1 ) V, 
\ "(n— m)x(m— 1) / 

where U, V are orthogonal matrices of appropriate dimensions. 

On the other hand, by Lemma^any pair ({/', V) of orthogonal matrices of appropriate size defines 
a doubly stochastic positive map 

M (f /',F')=(Í vúy,). 

Therefore 

P max = max F(M(U', V')) = max tr(V T D T U T P 2 U DV PA 

U',V U',V 

= maxtr(D T (U T P 2 U)D(VP 1 V T )). 

The pair (U, V) of orthogonal matrices maximizes the function F(M(U', V')). 

Denote VP\V T by Pi and U T P 2 U by P 2 . Then the first order maximality condition is given by 
the commutation relations [D T P 2 D, Pi] = \DP\D T , P 2 ] = 0. Partition the matrix P 2 as 

where P 2 n is of size (m — 1) x (m — 1). Then above commutation relations imply [Pi, P 2 1:L ] = 0, P 2 12 = 
Pi 1 = 0. Let now W m -\ be an orthogonal matrix that simultaneously block-diagonalizes Pi and P^ 1 , 
and let W n - m be an orthogonal matrix that diagonalizes P 2 2 . Now note that P max = tr(D T P 2 DP{) = 
ír(í> T [(íioj(W m _i, W n - rn ) P 2 diag(W m -i, W n - m ) T ]D[W m -iPi W^_ 1 \). The produets in brackets are 
diagonal and have the form (U') T P 2 U' , V Pi(V') T for some orthogonal matrices U', V. Hence we can 
assume without loss of generality that Pi , P 2 are both diagonal. Therefore there exist pairwise distinct 
indices j\, . . . ,j m -x <= {1> . . . ,n — 1} such that P max = J2T=i Afc(Pi)Aj fe (P 2 ). Obviously this sum is 
maximal if jk = k for all k, and we arrive at the inequality P max < X^l-Li 1 Afc(Pi)Afc(P 2 ). 

On the other hand, there exist orthogonal matrices U', V such that 
(U') T P 2 U' = diag(X 1 (P 2 ),X 2 (P 2 ), . . . , A„_i(P 2 )), VP^VY = diag{X x {Pxl . . . , A ro _i(Pi)). Then we 
have F(M(U',V')) = Y%=1 Afc(Pi)A fc (P 2 ) and P max > A)t(Pi)A fe (P 2 ). The proof is complete. 

□ 

We have proven the following 
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Corollary 4 Let r be the largest number such that the inclusion K^aii ( r ) C K sep holds. Then 



min(n,m) — 1 

max<j -l + (l + A 1 (P 1 ))(l + A 1 (P 2 )), A fc (Px)A fc (P 2 ) 



fc=i 



-1/2 



□ 



By Lemma[H]we now have the following theorem. 
Theorem 2 The radius of the largest K\ ® K^-separable ball around e™ ® eg is given by 

-1/2 



min(n,m) — 1 

max l (1 + A 1 (P 1 ))(1 + Ai(P 2 )), 1 + ^ A fc (Pi)A fc (P 2 ) 



fc=i 



□ 



Corollary 5 Let B\ C R m , P 2 C R" be balls of radií p\ : p2 < 1 around the unit vectors e™,eç, 
respectiu ely. Let K\, K2 be the cones generated by these balls. Then the radius of the largest K\ ® K2- 
separable ball around the unit vector e™ ® eg € R mn equals 



[max { Pl 2 p 2 2 , 1 + (min(n, m) ~ l){ Pl 2 - l)(p 2 2 - 1)}] 



-1/2 



The corollary is a direct consequence of the preceding theorem and Lemma|Hl 

5 Application to multi-qubit systems 

In this section we apply the obtained results to compute largest Ki g) X 2 -separable balls of bipartite 
matrices around the identity, where the cones K\ , K 2 are generated by balls around the identities in 
the factor spaces. We provide the exact value of the radius of such largest balls in dependence on the 
radii of the original balls and the dimensions of the matrices. These results will be used to compute 
lower bounds on the largest separable ball of unnormalized mixed states for multi-qubit systems. 

Denote the space of k x k hermitian matrices by H(k). Let B ri C 7í(m), B r2 C Ti.(n) be balls of 
radii n < \/m, r 2 < \fn around the corresponding identities J m , /„ and let K\, K% be the conic hulls 
of these balls. We look for the largest ball around the identity I nm 6 Ti{mn) = H(m) <g> Ti(n) which is 
contained in the cone of K\ ® ÜT 2 -separable matrices. 

The following corollary is a consequence of Corollary [S] and the fact that the identity in H(n) = R" 
has norm y/n. 

Corollary 6 The largest ball around L nm e TL(mn) = 7í(m) (£> TL{n) which is contained in the cone of 
K\ ® K2-separable matrices has radius 



r = mm ?'ir 2 , 



Imnr\r2 



y (min(m 2 , n 2 ) — l)(m — r\)(n — rf) + v\r. 



□ 



We see that for large dimensions and small ri,r 2 r is asymptotically equal to mi „^ n ^ ■ This 

asymptotics was independently found by Leonid Gurvits 1 . 

Let us use this result to obtain a bound on the radius of the largest separable ball of unnormalized 
density matrices for multi-qubit systems. Let m = 2, r% = 1 and set n{k) = 2 k ~ 1 . Define a sequence 
Pk recursively by p\ = 1 and 



■ . V 'mn{k)r x p k -x . 
p k = mm | npk-i, v = | (11) 

\J (min(m 2 , n(kf) - l)(m - r\)(n{k) - p^_ x ) + r 2 p 2 k _ K 
V2^pk-i 1 \/2V fe _i 



mm pk-i, 



1 Leonid Gurvits, personal communication 
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for k > 2. It follows that 




— 2 ^ — 2 

Pk = 2^-1 



and we get the explícit expression 



1(1 

3 ^2 



) 



k 



2 k/2 



Pk 



-2 



+ 2 



fe 



V3*- 1 + 1 



Theorem 3 pk 



2 k/2 



is a lower bound on the radius of the largest separable ball of unnormalized 



a/3 ! =- 1 +i 



multi-partite mixed states of a k-qubit system around the identity matrix in the space Tt(2)® k . 
Proof. We prové the theorem by induction. 

For a one-qubit system p\ = 1 is the radius of the largest ball around I2 in the cone 7í+(2) of 
positive scmidcfmitc hcrmitian 2x2 matrices. Hence for k — 1 the bound pk is exact. 

Assume now that the ball Bk-i C 7í(2)®( fc-1 ' of radius pu-i around the identity matrix I 2 h-i £ 
7Y(2)®' fe -i) consists of unnormalized separable states of a (fe — l)-qubit system. Let us apply Corollary[ü] 
with m = 2 and r± = 1. Since the cone 7í+(2) is isometric to the Standard Lorentz cone L4 = K st (I^), it 
will be generated by a ball of radius 1 around I2 and we get K\ = TÍ+(2). Let further n = n(k) — 2 k ~ 1 
and r2 — Pk-i- If we identify the space Tí(n) with the space 7í(2)®' ,c_1 \ then the cone K2 will be 
generated by Bk-i- 

But then the ball Bk C 7í(2)® fc of radius pk around the identity matrix J 2 fc is W+(2) ® Bk-\- 
separable by and Corollary [U Thus it is also (2)® -separable by the assumption on Bk-i- 
□ 

Remark: pk is the best bound one can obtain by tensoring in the spaces H(2) successively and 
approximating each time the separable cone by the largest ball-generated cone contained therein. This 
general approach was proposed and successfully applied by Gurvits and Barnum in j^j. 

Remark: Since both factor cones in the 7i+(2) l82 -separable cone are isometric to L4, CorollarylHl 
provides the exact result also for k — 2. Gurvits and Barnum obtained the exact result for a general 
bipartite space in 0]. 

For a 3-qubit system we get a radius of ^/4/5 instead of -y/8/11 and for n-qubit systems with n > 4 
an improvement of over 12.3% with respect to Gurvits' result in [Hj. The new bounds imply that with 
Standard NMR preparation technique one needs at least 36 qubits to obtain entanglement, which is a 
slightly stronger restriction than the one proven by Gurvits and Barnum 6 . 

6 Conclusion 

In this contribution we dealt with cones consisiting of elements separable with respect to two Lorentz 
cones. Such cones are prospective candidates for the approximation of more complex separable cones, 
such as the cones of unnormalized separable states of a multi-partite quantum system. The idea of 
using a Lorentz cone to approximate one of the factor cones in a bipartite setting and recursively 
in a multi-partite setting was introduced by Leonid Gurvits and Howard Barnum in [S]. Later they 
obtained asymptotically exact results on the size of largest separable balls in [ü] . 

We considered different aspects of ball-ball separable cones. Theorem Q]describes the extreme rays 
generating the cone dual to a ball-ball separable cone, i.e. a cone of Lorentz-to-Lorentz positive maps. 
There are two kinds of such rays, and all rays of one kind are equivalent under the action induced 
by the automorphism groups of the individual Lorentz cones. Correspondingly, the ball-ball separable 
cones possess two kinds of largest faces. Here by a "largest" face we mean a non-trivial face that is not 
the intersection of other, strictly larger faces. The shape of these faces is described in Propositions^ 
and El In CorollariesEl and|21 we established that the largest faces are highly intersecting each other, 
unlike the largest faces of a single Lorentz cone. In Theorem0and Corollary[S]we compute the radius 
of the largest ball around an element on the central ray of a ball-ball separable cone that is contained in 
this cone. This result extends to the case of balls in ellipsoid-ellipsoid separable cones. Such cones are 
affinely isomorphic, but not isometric to a ball-ball separable cone. The extension to ellipsoid-ellipsoid 
separable cones allows to use more flexible approximations of individual factor cones by ellipsoidal 
cones, which in may be more appropriate than Lorentz cones in some situations. 

Finally, we applied the developed theory to the case of a multi-qubit quantum system. Due to the 
exactness of our estimates we were able to sharpen previously available bounds on the radii of maximal 
separable balls around the uniformly mixed state. Our bounds in Theorem |21 are about 12% tighter 
than the best bounds obtained so far 0. 
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